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COMPUTATIONAL  FLUID  DYNAMICS  IN  THE  UNITED  KINGDOM* 

by 

M.  G.  Hall 
S.  P.  Fiddes 


SUMMARY 


A  review  is  presented  of  computational  methods  in  aerodynamic  research 
and  design,  with  application  of  the  methods  and  associated  computing  facilities 
included.  The  review  begins  with  a  brief  survey  of  the  field,  to  give  an  overall 
view  and  to  identify  what  seem  to  be  notable  features.  These  are  then  described 
in  turn.  The  first  is  a  development  by  Hall  and  his  colleagues  of  accurate  and 
fast  schemes  for  solving  the  Euler  equations,  based  on  the  finite-volume  cell- 
vertex  methods  introduced  by  Denton  and  Ni.  Next  is  a  development  of  a  block- 
structured  ('multiblock' )  grid  generation  technique  by  the  Aircraft  Research 
Association.  There  follows  an  application  of  this  grid  generation  scheme,  in 
conjunction  with  an  Euler  solver,  to  the  calculation  of  the  transonic  flow  past 
a  wing-body-canard  configuration.  Finally,  an  application  of  viscous-inviscid 
interaction  techniques  to  the  design  of  a  fan  rotor  is  described. 


*  Invited  AIAA  Paper  No.  87—1 1 32— CP ,  presented  at  AIAA  8th  Computational  Fluid 
Dynamics  Conference,  Honolulu,  USA,  9-11  June  1987. 
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Introduction 

The  review  begins  with  a  survey  that  13 
restricted  to  the  range  of  computational  methods 
that  are  in  use  or  under  active  development  for 
application  In  aerodynamic  research  and  design. 

The  survey  deals  first  with  methods  thaC  are  well 
established  but  for  which  no  major  developments 
are  envisaged,  and  then  moves  on  to  methods  that 
are  In  use  and  still  undergoing  active  develop¬ 
ment.  The  latter  fall  mainly  Into  four  classes  - 
methods  for  the  Euler  equations,  for  viscous- 
lnviscld  Interactions,  for  the  Navler-Stokes 
equations  and  for  grid  generation  -  and  these  are 
covered  In  turn.  For  each  class  of  method  some 
relevant  applications  are  noced.  In  the  course  of 
the  survey  a  number  of  features  are  singled  out 
for  more  detailed  description  in  Che  remainder  of 
the  review.  The  survey  concludes  with  a  brief 
accounc  of  computing  facilities. 


(TSP)  equacion  are  still  in  use.  However  Che 
unreliability  of  TSP  methods  when  perturbations 
are  not  small,  and  of  all  potential  methods  when 
shock  waves  are  present,  no  longer  oucweighs  chelr 
relative  simplicity.  Ic  seems  appropriate,  as  our 
TSP  methods  are  being  retired,  to  note  the  key 
role  they  have  played  in  the  history  of  CFD. 

Their  Introduction  by  Murman  and  Cole  and  their 
revolutionary  effect  on  Che  development  of  aero¬ 
dynamics  were  discussed  in  an  earlier  review  given 
at  the  5th  AIAA  CFD  Conference  in  19814  .  Since 
that  time  a  number  of  aircraft  designed  with  the 
aid  of  such  methods  have  gone  into  service.  For 
Che  UK  the  first  of  the  'CFD  airliners'  was  the 
A310  Airbus.  The  Impact  on  the  A310  of  the  avail¬ 
ability  of  CFD  methods  chat  could  adequately  cal¬ 
culate  flows  with  shock  waves  about  wing-body 
combinations  has  been  summarised  by  Jupp5 .  He 
notes  that  in  comparison  with  the  Inferior 
('pre-CFD')  A300  wing  design,  the  A3I0  design 
required  one-third  the  number  of  wlndtunnel  models 
and  13  months  less  time  to  arrive  at  the  final 
design.  The  Importance  of  CFD  in  the  design  of 
the  A320  wing  Is  similarly  described  by  Back  and 
Wedderspoon6 .  The  A320  Is  proving  an  outstanding 
commercial  success,  with  over  400  orders  placed 
before  the  first  flight. 

Another,  perhaps  less  well  known,  field 
where  potential  methods  have  had  a  revolutionary 
role  Is  that  of  propeller  design.  We  note  here 
the  radical  Improvements  that  have  been  obtained 
following  work  at  the  Aircraft  Research 
Association  (ARA)  on  the  design  of  high- 
performance  propellers  for  commuter  aircraft7 >8. 

There  are,  however,  certain  special  roles  in 
which  potential  methods  seem  likely  to  remain  use¬ 
ful  for  some  time.  Where  trends  rather  chan 
absolute  accuracy  are  important  (such  as  in  pre¬ 
liminary  design  studies,  or  in  unsteady  calcu¬ 
lations  for  flutter  estimation)  the  simple  TSP 
mechods  will  remain  useful.  On  the  other  hand, 


In  spite  of  the  restriction  to  aerodynamics 
the  range  of  methods  covered  here  is  so  wide  thac 
no  method  can  be  discussed  In  depth.  Reference 
may  be  made  to  a  number  of  specialised  reviews  by 
8rltlsh  auchocs  that  give  a  more  comprehensive 
accounc  of  cheir  subject  maccer.  In  particular 
the  reader's  attention  is  drawn  to  the  recent 
review  by  Dunham1  on  turbomachinery  applications, 
the  decailed  discussion  of  vlscous-inviscld  inter¬ 
actions  by  Lock  and  Williams2  ,  and  the  survey  by 
Roe3  of  characcerlstic-based  schemes  for  solution 
of  che  Euler  equations. 

Survey  of  Field 
Established  Mechods 

Methods  for  solving  Che  f ul 1-potent la  1 
e.uaclon  and  che  transonic  sraa Ll-perturbatlon 


f ull-potencial  mechods  offer  high  accuracy  for 
subcritlcal  flows  and  thus  can  provide  datum 
solutions  against  which  the  newer  Euler  mechods 
can  be  judged.  Panel  mechods  will  remain  useful 
in  the  calculation  of  low-speed  flows,  while 
field-panel  mechods  continue  to  offer  relatively 
simple  solutions  for  mildly  supercritical  flows 
past  complex  configurations.  Finally,  potential 
mechods  will  continue  to  have  a  role  in  composite 
schemes,  such  as  zonal  mechods  and  vlscous- 
inviscld  interaction  schemes,  where  che  potential 
equacion  is  valid  for  an  outer  part  of  che  flow. 

Another  type  of  mechod  chat  provides  a  high 
standard  of  solution  is  che  space-marching  scheme 
for  the  Euler  equations.  The  limitation,  wich 
such  schemes,  Co  supersonic  flow  everywhere  is 
leading  to  che  predominance  of  time-marching  Euler 
mechods  for  practical  applications  in  the  UK.  At 
present,  however,  it  is  difficult  to  assess  the 
accuracy  of  the  latter  whereas,  for  example, 
accurate  solutions  have  been  obtained  for  the  wave 
drag  of  a  delca  wing  calculated  via  a  space¬ 
marching  method3. 
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Methods  for  Che  Euler  Equations 

Jameson's  original  cell-centred  finite- 
volume  scheme  for  the  Euler  equations  has  been 
developed  by  British  Aerospace  (BAe)  for  extensive 
practical  applications.  The  work  involved,  and 
some  applications,  have  been  described  by  Doe, 
Brown  and  Pagano10,  while  ocher  UK  applications, 
for  military  aircraft,  have  been  described  by 
Marchbank11.  An  important  development  of  the  work 
has  been  its  integration  with  the  ASA  multiblock 
grid  generation  technique. 

Another  feacure  of  UK  work  is  the  develop¬ 
ment  of  mulcigrld  finite-volume  schemes  of  the 
cell-vertex  type.  These  techniques,  as  originally 
proposed  by  Denton  and  Ni  suffered  from  short¬ 
comings  in  accuracy  and  speed  (see,  for  example, 
references  12  and  13),  but  Hall14  and  his 
colleagues  at  RAE  have  achieved  substantial 
Improvements  in  these  two  key  areas  and  have  pro¬ 
ceeded  on  to  extensions  of  Che  basic  multigrid 
method15,  che  inclusion  of  adaptive  grids1®,  and 
the  development  of  a  three-dimensional  method17. 
Morton  and  Paisley18  and  Roe19  have  placed  such 
schemes  on  a  more  secure  mathematical  foundation. 
Reference  18  also  includes  results  for  two- 
dimensional  cell-vertex  Euler  solutions  with 
fitted  shocks.  An  account  of  these  develop¬ 
ments  of  che  cell-vertex  mechod  will  be  presented 
below. 

Vlscous-Inviscld  Interactions 

UK  work  on  viscous-lnviscid  interactions  is 
described  in  the  detailed  and  wide-ranging  review 
by  Lock  and  Williams2.  The  UK  continues  to  favour 
integral  methods  for  boundary  layer  calculations 
mainly  because  of  their  speed  and  simplicity. 

They  have  also  been  developed  to  yield 
sufficiently  accurace  solutions  for  a  useful  range 
of  practical  problems.  Developments  in  techniques 
for  coupling  boundary  layer  methods  to  inviscid 
methods  (both  Euler  and  potential)  have  produced 
methods  chat  can  predict  flows  involving  boundary 
layer  separation,  with  most  success  being  achieved 
for  two-dimensional  flows  thus  far.  In  three 
dimensions  simple,  direct  coupling  has  been  imple¬ 
mented  in  a  number  of  mechods,  while  for  two- 
dimensional  methods  direct,  inverse,  semi-inverse 
and  qua3l-slmulcaneous  coupling  procedures  have 
all  been  used.  These  types  of  coupling  and  their 
relative  merits  are  discussed  In  Che  above  review. 
Cross29  gives  a  description  of  a  quasi- 
slmulcaneous  technique  and  its  use  in  conjunction 
with  a  panel  mechod.  An  example  of  the  applica¬ 
tion  of  coupled  vtscous-invlscid  methods  to  the 
design  of  a  fan  rotor  13  described  by  Cinder  and 
Calverc21,  and  chls  is  discussed  below  as  one  of 
our  chosen  highlighcs. 

Methods  for  Navler-S  Cokes  Equations 

Although  there  Is  considerable  research 
effort  in  the  UK  relating  to  the  Navier-Stokes 
equations,  there  are  relatively  few  methods  that 
are  in  use  or  being  actively  developed  for  aero¬ 
dynamic  research  and  design.  The  mechods  that 
have  found  application  are  based  on  mechods 
developed  outside  che  UK.  Thus  Rolls-Royce  have 
under  continuing  development  three-dimensional 
methods  based  on  the  work  of  J.  and  J.G.  Moore22 
and  Thorapklns  and  his  colleagues21  In  USA.  These 
methods  have  been  used  as  analysis  tools,  in  Che 


design  of  che  high-pressure  curblne  nozzle  guide 
vanes  for  che  RB211-535  and  RB211-524,  and  in  che 
design  of  the  compressor  ouclec  guide  vanes  for 
che  IAE-V25G0.  Progress  in  che  development  of 
practical  Navier-Stokes  mechods  is  presently 
hampered  not  only  by  inadequacies  in  computing 
capacity,  but  also  by  fundamental  difficulties: 

Che  generation  of  suitable  grids,  the  modelling  of 
turbulence,  and  the  attainment  of  satisfactory 
accuracy  and  convergence  with  grids  thac  are  very 
highly  stretched  or  distorted. 

Grid  Generation 

Grid  generation  methods  for  complex  con¬ 
figurations  make  up  che  last  class  of  methods  that 
are  both  in  use  and  under  acclve  development. 
Because  of  their  practical  Importance,  several 
different  approaches  have  been  explored.  These 
include  non-aligned  grids24'25,  where  the  grid  is 
not  aligned  with  the  body  surface  (leading  to 
simple  grid  generation  but  complicated  treatment 
of  boundary  conditions),  unstructured  grids,  over¬ 
lapping  grids,  embedded  grids28,  solution-adaptive 
grids1®  and  block  structured,  or  multiblock, 
grids27*28.  The  power  of  the  multiblock  approach 
has  encouraged  UK  researchers  to  puc  considerable 
effort  into  the  development,  evaluation  and  prac¬ 
tical  use  of  multiblock  mechods.  A  multlblock 
framework  can  incorporate  the  ocher  types  of  grid 
and  be  applied  to  most  (if  not  all)  practical 
configurations.  The  computational  advantages  of 
having  a  regular  structured  grid  can  be  retained 
for  individual  blocks,  although  the  array  of 
blocks  themselves  will  in  general  be  unstructured. 
The  development  of  Che  method  has  involved  study 
of  several  questions:  how  to  define  individual 
blocks  and  to  set  up  che  array  of  blocks  to  cover 
the  complete  field;  what  degree  of  continuity  of 
grid  lines  at  block  interfaces  is  optimal;  how  to 
control  the  point  or  cell  density;  how  to  avoid 
singularities;  and  what  degree  of  automation  of 
grid  generation  is  feasible,  or  even  desirable? 

An  account  will  be  given  here  of  the  multiblock 
mechod  developed  at  the  Aircraft  Research 
Association  (ARA),  and  its  application  to  com¬ 
puting  Che  flow  past  a  wing-body-canard 
configuration. 

Computing  Facilities 

The  computing  facilities  used  for  CFD  in  the 
UK  range  from  micro  and  mini-compucers  and 
graphics  workstations  through  conventional  main¬ 
frames  to  supercomputers.  The  main  users  of 
supercomputers  are  groups  in  the  airframe  and 
aero-engine  industry,  who  access  such  facilities 
through  bureaus.  Other  major  supercomputer  users 
are  RAE  and  ARA  (using  the  Cray  1-S  sited  ac  RAE 
Farnborough)  and  Universities  (who  have  network 
access  to  a  Cyber  205,  a  Cray  X-MP  and  2  Cray 
l's).  The  larger  institutions  have  their  own 
mainframe  computers,  but  the  trend  is  towards 
greater  use  of  supercomputers  for  CFD.  Associated 
with  this  trend  is  a  rapid  growth  in  the 
facilities  available  for  pre-  and  post-processing, 
particularly  graphics  workstations.  There  is  a 
growing  awareness  of  Che  need  for  more  fast  memory 
in  supercomputers.  This  need  may  not  be  met  by 
out-of-core  memory  unless  it  is  both  efficient  and 
transparent  in  tcs  use. 

The  latesc  supercomputers  are  characterised 
by  having  a  small  number  (typically  four)  of  very 
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high  speed  vector  processors  sharing  a  large 
amounc  (megawords)  of  common  memory.  These  offer 
limited  prospects  for  multiprocessing,  boch  In  the 
simpler  single-inscrucclon  multiple-data  (SIMD) 
mode  where  the  processors  execuce  the  same 
instruction  simultaneously  on  different  data,  and 
in  the  multiple-instruction  (MIMD)  mode  where, 
even  with  only  a  few  processors,  the  programming 
is  far  from  trivial.  However,  in  the  United 
Kingdom,  the  concept  of  aultlprocessing  hardware 
has  been  taken  much  further.  Tw<^  noteworthy 
developments  will  be  discussed  -  the  Distributed 
Array  Processor  (DAP)  originally  developed  by 
International  Computers  Limited  and  the  Transputer 
developed  by  Inmos. 

The  DAP  is  a  SIMD  device  intended  for  use 
wlch  a  host  computer.  By  analogy  with  the  concept 
of  a  veccor  processor,  the  DAP  may  be  viewed  as  a 
matrix  processor.  The  processor  array  currently 
comprises  1024  single-bit  processors  arranged  as  a 
32  by  32  array.  Each  processor  is  supported  by  a 
32K  bit  memory.  Communication  between  processors 
is  local,  to  the  four  nearest  neighbours.  Grosch28 
has  recently  described  the  use  of  a  DAP  in  con¬ 
junction  wich  a  two-dimensional  Navier-Stokes 
solution  method  based  on  relaxation.  The  archi¬ 
tecture  of  the  DAP,  as  a  matrix  processor,  is  well 
suited  to  such  an  approach.  However,  the  present 
DAP  architecture  is  rather  limited  in  size  and 
processing  power,  and  substantial  increases  in 
memory,  speed  and  communications  capability  of  a 
DAP-like  device  are  required  before  practical 
three-dimensional  problems  can  be  contemplated. 


which  would  be  used  all  of  the  time.  A  more 
realistic  approach  may  be  to  consider  'super- 
blocks',  where  a  group  of  blocks  are  associated 
wich  each  processor.  In  any  case  mulclblock  seems 
to  offer  a  nacural  ' task-splicclng '  scheme  for  use 
of  a  multiprocessor  syscem  in  the  calculation  of 
flowflelds  past  complex  configurations. 

Cell-Vertex  Multigrid  Solution  of  the 
Euler  Equations 

Cell-Vertex  Formulation 

The  cell-vertex  formulation14  differs  from 
the  better  known  cell-centred  formulation  in 
having  the  dependent  variables  specified  ac  the 
vertices  of  a  computational  cell  rather  than  at 
the  centre.  In  both  formulations  the  fluxes 
through  the  sides  of  the  cell  are  approximated  in 
order  to  evaluate  the  flux  balances  and  hence  the 
changes  in  dependent  variables  for  the  entire 
cell.  The  approximation  used  in  the  cell-vertex 
scheme  for  the  flux  through  a  cell  side  is  the 
trapezoidal  rule  for  integration  and  hence  is 
second-order  accurate  in  cell  size  whacever  che 
grid,  whereas  the  corresponding  cell-cencred 
approximation,  being  an  average  of  quantities  ac 
the  centres  of  adjacent  cells,  will  only  be 
second-order  accurate  on  nearly  uniform  grids.  An 
additional  advantage  of  the  cell-vertex 
formulation  is  that  no  extrapolation  is  required 
to  obtain  values  for  the  dependent  variables  on 
body  surfaces. 


By  contrast,  the  Transputer  (which  is  to  be 
regarded  as  an  element  of  a  larger  system)  com¬ 
prises  a  powerful  32-bit  processor  giving  10  MIPS 
processing  speed,  2  kilobyte  of  memory  and  four 
links  each  providing  a  10  megabit  per  second  data 
highway,  ie  it  does  not  use  a  bus  structure  for 
data  transfer  between  processors.  It  comes  on  a 
single  chip  and  is  a  useful  computer  in  its  own 
right.  However,  ic  Is  designed  to  be  a  building 
block  for  a  true  multiprocessor  system  in  that  a 
Transputer  system  can  be  configured  to  the  user's 
requirements  and  operaced  in  either  SIMD  or  MIMD 
mode,  or  in  both  modes  simultaneously.  Recent 
developments  Include  che  ability  to  reconfigure 
the  syscem  by  software. 

These  impressive  developments  in  hardware 
and  computing  architecture  now  pose  a  challenge 
for  che  CFD  mechod  developer  -  how  to  use  them. 

The  burden  of  using  a  multiprocessor  system 
effectively  rest3  with  the  applications  program¬ 
mer  -  there  are,  as  yet,  virtually  no  software 
tools  chat  will  cake  a  code  developed  on  a 
'general  purpose’  computer  and  partition  che  code 
for  execution  over  a  number  of  processors  -  the 
programmer  muse  perform  this  'task  splitting' 
himself.  Ac  presenc  chls  is  not  a  desirable  or 
practical  computing  environment  for  code 
developers  concerned  with  developing  new  methods 
or  crying  new  applications.  However,  it  may  be 
feaslbLe  co  implement  on  a  multiprocessor  syscem  a 
proven  code  chat  is  used  on  a  production  basis. 

The  multiblock  concepc  suggests  a  possible  way 
forward  in  this,  where  each  processor  could  calcu¬ 
late  che  flow  in  one  block.  Note,  however,  that 
the  wing-body-cana rd  configuration  referred  co 
above  uses  430  blocks,  and  a  'general-purpose' 
multiblock  multiprocessor  syscem  may  therefore 
require  several  hundred  processors,  not  all  of 


Recently,  studies  of  accuracy  by  Morcon  and 
Paisley18  and  Roe19  have  yielded  Che  following 
results.  For  a  quadrilateral  cell  with  sides 
0(h)  and  flux  balance  0(h2)  the  truncation 
error  in  the  flux  balance  is  0(h3 )  or  first 
order  in  the  cell-vertex  formulation.  Irrespective 
of  the  size  of  the  adjacent  cells.  For  a  quadri¬ 
lateral  with  opposite  sides  differing  by  Q(h2  )  , 
approximating  a  parallelogram,  the  truncation 
error  is  0(h4  )  or  second  order.  Such  quadri¬ 
laterals  are  frequently  found  in  practice,  for 
example  when  che  grid  is  generated  analytically, 
or  by  repeated  bisection  of  an  arbitrary  coarse 
grid.  In  contrast,  the  truncation  error  of  che 
cell  flux  balance  on  Che  cell-centred  formulation 
depends  on  the  sizes  and  shapes  of  che  adjoining 
cells.  Analysis  in  this  case  is  difficult.  How¬ 
ever,  for  a  non-uniform  reccangular  grid  (for 
which  che  cell-vertex  flux  balance  is  second-order 
accurate)  the  racio  of  successive  cell  sides  must 
be  [ (1  +  0(h)] :1  for  first-order  accuracy  and 
[ (1  +  0( h2  )  ] : 1  for  second-order  accuracy. 
Furthermore,  if  che  vertex  of  a  cell  of  an 
initially  rectangular  grid  is  displaced  by  an 
amount  e  ,  Chen  E  must  be  0(h2  )  for  first-order 
accuracy,  while  e  of  order  h3  is  required  for 
second-order  accuracy.  (By  contrast,  the  cell- 
vertex  scheme  yields  second-order  accuracy  for  E 
of  order  h2  .)  The  truncation  error  arising  from 
Che  simple  averaging  of  quantities  at  the  centre 
of  adjacent  cells  can  be  reduced  (in  principle)  by 
weighting  the  average.  However  this  complicates 
che  formulation  and,  more  seriously,  makes  ic  more 
difficult  to  ensure  the  convergence  of  the 
iterative  process- 

Spurious  modes  are  anocher  possible  source 
of  error.  Analysis18  shows  that  for  a  model  two- 
dimensional  problem  the  cell-vertex  and  cell- 
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centred  formulations  yield  the  following  spurious 
modes: 

+  -  + 

Cell-vertex:  -  +•  - 

+  -  + 

+  -+  +  -  +  +  +  ■*• 

Cell-centred:  -  +  -  +-+  —  —  - 

+  -  +  +  -  +  +  +  + 

Whilst  the  cell-vertex  scheme  has  one 
spurious  mode  (the  well-known  chequer-board  mode), 
the  cell-centred  scheme  has  three  possible  modes. 
Thus  it  should  be  expected  that  the  prevention  and 
damping  of  spurious  modes  will  be  easier  in  the 
cell-vertex  than  in  the  cell-centred  formulation. 

Iteration  to  Steady  State 

In  the  original  version  of  Hall's  cell- 
vertex  multigrid  method  the  solution  was  Iterated 
to  the  steady  state  by  use  of  a  Lax-Wendroff  time- 
stepping  algorithm  in  combination  with  a  simple 
multigrid  algorithm.  The  emphasis  during  the 
development  of  this  early  method  had  been  feasibi¬ 
lity  and  accuracy,  rather  than  speed.  Two  devices 
for  increasing  speed  have  since  been  incorporated, 
and  are  now  described. 

Consider  first  the  changes  calculated  on  the 
coarser  grids  of  a  multigrid  cycle.  In  the  origi¬ 
nal  method  the  changes  were  simply  a  redistribu¬ 
tion  and  re-scaling  of  changes  calculated  on  Che 
previous  finer  grid.  It  would  be  expected  that  a 
direct  integration  of  che  Euler  equations  (in  a 
form  appropriate  for  defining  changes  in  the  solu¬ 
tion)  would  accelerate  the  rate  of  convergence  to 
a  steady  stace  -  indeed,  Jameson  adopts  such  a 
technique.  Thus  a  time-stepping  algorithm  of  Che 
Runge-Kucca  type  was  introduced'-'*  ,  for  calcula¬ 
tions  on  the  coarser  grids  only.  It  was  tried  at 
all  grid  levels  but  found  then  to  be  slower.  The 
algorithm  differs  from  Jameson's,  however,  in  that 
lc  does  not  include  fourth-order  spatial  dissipa¬ 
tion,  but  includes  a  temporal  damping  provided  by 
part  of  the  Lax-Wendroff  increment.  This  adds  no 
additional  cruncacton  error  co  che  sceady-scace 
solution  and  is  also  more  compact. 

The  second  accelerative  device  i3  a  simple 
form  of  residual  averaging,  applied  on  the  coarser 
grids  only.  In  che  original  method  changes  on  the 
coarser  grids  were  interpolated  linearly  to  give 
the  changes  at  each  of  the  points  on  the  finest 
grid.  The  interpolation  is  now  smoothed,  thus 
averaging  che  residuals. 

The  convergence  and  stability  of  the  above 
scheme  have  been  scudied  by  Morcon  and  Paisley'®. 
With  the  cell-vertex  formulation  and  Lax-Wendroff 
clrae-scepplng  it  is  shown  that  Ni's  simple,  non- 
reflectlve  boundary  condition  (in  which  any  resi¬ 
duals  lying  outside  the  computational  domain  are 
set  co  zero)  is  a  necessary  condition  for  con¬ 
vergence.  In  cell-centred  formulations  there  is 
no  simple  boundary  treatment.  Morton  and  Paisley 
deduce  thac  provided  che  boundaries  are  treated 
correctly  che  residuals  in  the  individual  cells 
will  be  decoupled  so  chat  at  convergence  all  che 
individual  residuals  will  vanish,  rather  than  any 
averages.  Non-uniformity  of  a  grid  is  shown  co  be 
destabilising.  Hall  had  found  it  effective,  in 
averaging  che  residuals  from  four  neighbouring 


cells  to  obcaln  an  updace  for  che  common  vertex, 
to  weighc  che  residuals  by  individual  cell  areas. 
This  weighted  average  is  shown  to  emerge  naturally 
from  an  alternative  f lnite-elemenc  formulation  of 
Che  Taylor-Galerkin  cype.  Multi-stage  Runge-Kucca 
algorithms  are  shown  co  be  destabilised  by  non- 
uniformity.  With  che  cell-vertex  formulation  this 
can  be  alleviated  by  weighting  residuals,  but 
further  damping  may  be  required.  With  che  cell- 
centred  formulation  the  addition  of  excra  dis¬ 
sipation  can  stabilise  the  iteration,  and  Indeed 
the  foureh-order  dissipation  of  Jameson  is  a 
popular  remedy. 

Shock  Fitting 

Although  the  cell-vercex  formulation  has 
important  advantages  in  accuracy  and  in  Che  number 
of  spurious  modes,  it  will  still  be  subject  to 
significant  errors  if  shock  waves  are  captured,  in 
common  with  other  formulations.  Morton  and 
Paisley18  have  developed  a  shock-fitting  scheme 
and  thus  allow  the  cell-vertex  scheme  to  yield 
accurate  solutions  for  supercritical  flows  with 
embedded  shock  waves.  The  scheme  involves  an 
adaptive  shock-aligned  grid.  The  conventional 
shock-capturing  version  of  the  cell-vertex  scheme 
is  used  to  give  an  initial  estimate  of  the  shock 
position.  The  grid  is  then  aligned  with  the 
estimated  shock  locacion,  and  che  grid  line 
running  along  the  shock  is  treated  as  a  boundary 
of  the  computational  domain.  A  cell  on  the 
upstream  side  of  the  shock  is  treated  as  a  super¬ 
sonic  outflow  boundary  cell,  while  a  cell  on  the 
downstream  side  is  treated  as  a  subsonic  inflow 
cell,  with  the  Rankine-Hugoniot  conditions 
replacing  the  usual  inflow  matching  conditions. 

At  each  time  step  the  calculation  yields  a  shock 
speed.  The  shock  is  moved  accordingly  and  the 
grid  re-aligned  repeatedly  until  a  steady  state  is 
reached. 


Results  for  NACA  0012  Aerofoil 

A  selection  of  results  for  Che  NACA  0012 
aerofoil  is  presented  Co  illuscrace  Che  accuracy 
and  speed  of  the  cell-vertex  multigrid  mechod 
described  above.  Figure  1  shows  Mach  number 
distributions  on  che  aerofoil  surface  for  captured 
and  fitted  shocks.  Table  1  shows  a  comparison  of 
calculated  lift  coefficients  for  cell-vertex  and 
cell-centred  oechods,  for  chree  cest  cases  on  boch 
coarse  and  fine  grids30.  The  same  C-mesh  was  used 

Table  1  Calculated  lift  coefficients.  NACA  0012 


Case 

No. 

cells  — 

Captured 

shocks 

Ceil  Cell 

centred  vertex 

Fitted 

shocks 

cell 

vertex 

Other 

results 

M»-0.30 

128»  16 

1.23 

1.271 

1.2733 1 

o-lO* 

256*32 

1.26 

1.273 

M.-0.80 

12Sx  16 

0.337 

0.366 

0.360 

0.36232 

a-1.25* 

256- 32 

0.365 

0.353 

M„-0.S5 

128x  16 

0.401 

0.380 

0.382 

0.39432 

o-l  .0" 

256«  32 

0.457 

0.406 
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for  boch  sets  of  shock-capturing  calculations. 
Included  la  the  table  are  results  obtained  by 
cell-vertex  shock-fitting  and  results  from  other 
accurate  calculations.  For  the  first  test  case  an 
accurate  potential  flow  solution  has  been  obtained 
by  Salas  et  al31;  it  can  be  seen  chat  the  corre- 
spondlag  cell-vertex  results  are  in  very  good 
agreement  with  che  potential  flow  result.  It  can 
also  be  seen  thac  the  results  from  use  of  shock 
fitting,  in  the  second  and  third  tesc  cases,  agree 
well  with  the  results  obtained  by  Pulliam  and 
Barton32  with  shock  capturing  on  an  exceedingly 
fine  grid.  The  results  from  the  cell-vertex 
method,  with  captured  shocks,  are  closer  to  the 
accucace  results  than  the  results  from  the  cell- 
centred  method.  Moreover,  che  cell-vertex  results 
obtained  on  the  coarse  grid  are  close  to  those 
from  the  fine  grid. 

Table  2  shows  a  comparison  of  CPU  times 
(Cray  1-S)  for  the  lift  to  reach  within  0.25?  of 
the  fully  converged  value.  The  benefits  of 
Including  Runge-Kutta  time-stepping  and  residual 
averaging  in  the  cell-vertex  method  are  clear. 

The  results  from  Jameson's  1985  cell-centred 
method33  (which  included  both  the  above  acceler¬ 
ative  devices)  were  obtained  on  the  same  grid  as 
used  for  the  cell-vertex  calculations. 


A  three-dimensional  version  of  the  cell- 
vertex  multigrid  method  for  swept  wings,  based  on 
the  SAE  work  for  two-dimensional  solutions,  has 
been  developed  by  Salraond17.  Only  Lax-Wendroff 
clme-scepplng  was  used  in  this  scheme.  More 
recently  che  method  has  been  accelerated  by  incor¬ 
porating  Runge-Kutta  time-stepping  and  smooched 
interpolation  (residual  averaging),  and  a  typical 
CPU  time  (Cray  1-S)  for  the  lift  to  converge  to 
0-52  of  che  final  value  is  140  seconds  for  a 
128  «  16  x  24  grid.  Anocher  Improvement  has  been 
che  incorporation  of  an  asymptodc  far-field 
boundary  condition  that  more  correctly  represents 
che  behaviour  of  a  lifting  solution  in  che  far 
field.  Klunker's31*  transonic  small-disturbance 
approximation  for  the  far-field  solution  is  used. 
This  involves  an  integral  over  the  surface  of  the 
wing  co  determine  che  appropriate  condition,  for 
every  point  on  the  far-field  boundary.  Table  3 
shows,  for  Che  ONERA  M6  wing,  Che  variation  of 
cocal  lift  with  distance  of  che  far-field  boundary 
from  the  wing.  It  can  be  seen  chat  macchlng  the 
boundary  conditions  to  a  small-disturbance  far- 
field  solution,  rather  than  free-stream  flow, 
allows  the  boundary  to  be  brought  slgnlf tcancly 
closer  co  che  wing  without  loss  of  accuracy. 


Table  3  Variation  of  lift  wich  boundary  dlscance. 
ONERA  M6  wing,  M.-0.84,  ct-3.06“ 


Distance 
chords  (spans) 

Free-stream 

condition 

Small-disturbance 

condition 

2.5 

(0.84) 

0.271 

0.287 

5 

(1.69) 

0.285 

0.290 

10 

(3.37) 

0.286 

0.288 
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Introduction 

The  guiding  principle  of  che  mulciblock  grid 
generation  method  developed  at  ARA27>2®  is  thac 
with  each  component  of  the  configuration  about 
which  a  grid  is  to  be  generated  there  is  associ¬ 
ated  a  particular  type  of  grid  thac  is  best  suited 
to  describing  the  flow  past  that  component.  This 
contrasts  with  the  more  common  approach  where  a 
block-based  structure  is  seen  primarily  as  a  means 
of  reducing  the  in-core  memory  requirements  of  a 
method.  The  ARA  multiblock  method  can  be  used  to 
generate  H-grids  everywhere  in  che  field,  but  It 
can  do  far  more;  in  particular,  it  offers  a  capa¬ 
bility  for  automated  generation  of  component 
adaptive  grids.  The  price  paid  for  such  flexi¬ 
bility  and  power  i3  complicated  program  logic. 

The  blocks  themselves  are  simple,  non-overlapping 
volumes,  each  with  six  faces  and  eight  corners.  A 
key  feature  is  Chat  only  one  type  of  boundary  con¬ 
dition  (eg  solid  boundary,  specified  transpiration 
velocity,  continuity,  etc)  is  permitted  on  each 
face  of  each  block.  This  simplifies  the  logic  but 
Increases  che  number  of  blocks  required  to  cover 
the  flow  field. 

The  grid  generation  procedure  can  be  broken 
down  into  three  stages:  topology  definition, 
surface  grid  generation  and  field  grid  generation. 
Successful  topology  definition  is  crucial.  A  good 
topology  can  ease  the  grid  control  problems  both 
on  the  surface  and  in  the  field.  An  inappropriate 
topology  can  make  It  impossible  to  generate  a 
satisfactory  grid.  Each  of  the  three  scages  are 
now  described  In  turn.  The  accounc  of  mulciblock 
concludes  with  a  discussion  of  the  prospects  for 
automated  grid  generation. 

Topology  Definition 

The  block  topology  is  the  arrangement  of  che 
blocks  with  respect  to  each  other.  Given  an  arbi¬ 
trary  configuration  the  definition  of  an  appro¬ 
priate  topology  is  a  challenging  task.  Inspection 
and  trial,  aided  by  experience  and  interactive 
computer  graphics,  is  an  obvious  strategy. 

Clearly,  a  systematic  approach,  with  access  co  a 
library  of  successful  topologies  and  at  least  some 
aucomacion,  is  desirable.  The  presenc  method 
requires  the  setting  up  of  a  'topology  file'  which 
lists,  for  each  block  face,  che  type  of  boundary 
condition,  the  adjacent  block,  che  adjacenc  face 
of  chat  block  and  che  oriencaclon  of  che  adjacenc 
block  with  respect  co  che  currenc  face.  This  file 
controls  the  complete  mulciblock  procedure. 

A  guide  to  che  construction  of  a  topology 
may  be  obtained  from  study  of  Figure  2.  Three 
basic  units,  "0",  "C"  and  ”H"  grids,  are  adopted 


Table  2  CPUsec  for  lift  to  converge  to  0.252 


Hall 

Hall 

Jameson 

No.  cells 

multigrid 

multigrid+ 

multlgrid+ 

1985 

1986 

1985 

128x16 

7.3 

6.3 

7.7 

256x32 

29.5 

13.3 

24.7 

Extension  Co  Three-Dimensional  Flows 
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for  the  construction  of  a  multibiock  grid.  O-grlds 
are  computationally  efficient,  but  present  topo¬ 
logical  difficulties  in  the  multiblock  framework. 
H-grids  are  the  least  efficient  computationally 
but  topologically  the  simplest.  C-grids  lie  in 
between  and  are  particularly  well  suited  for  flows 
past  wings.  The  figure  shows  two  possible 
topologies  for  a  simple  wing-body  combination: 

(i)  C-H-H,  ie  C-grid  around  the  wing  section, 
H-grid  spanwise  over  the  wing  tip,  and  H-grid 
vertically  around  the  body  section. 

(ii)  C-H-0 

Topology  (i)  is  simple  enough  to  be  constructed  by 
inspection.  The  section  AA'  shows  the  H-grid 
around  the  body  section,  and  the  (spanwise)  H-grid 
over  the  tip-  It  also  shows  how  the  wing  thick¬ 
ness  is  accommodated  by  opening  up  or  splitting  a 
surface  on  the  H-grid.  Sections  BB'  and  CC'  show 
how  the  wing,  and  its  extension  beyond  the  tip,  is 
embedded  in  a  C-grid,  and  how  the  C-grid  is  in 
turn  embedded  in  an  H-grid.  A  count  of  the  number 
of  blocks  in  the  BB'  and  CC'  sections  gives  the 
coral,  24,  for  the  complete  topology. 

Topology  (ii)  can  now  be  constructed.  The 
sections  CC'  and  DD'  are  identical  to  BB'  and  CC' 
respectively  in  topology  (i).  However,  the  intro¬ 
duction  of  an  O-grld  around  the  body  complicates 
the  topology  and,  as  shown  in  sections  AA'  and 
BB',  care  is  needed  in  interfacing  the  O-grid  with 
the  C-grid  over  the  wing  and  with  the  H-grid  above 
and  below  the  body.  Here  the  total  number  of 
blocks  is  42.  It  is  apparent  from  considering 
such  a  simple  configuration  that  the  policy  of 
admitting  only  one  type  of  boundary  condition  on 
each  block  face  does  indeed  lead  to  a  large  number 
of  blocks. 

Surface  Grid  Generation 

Suitable  grids  are  generated  on  the  surface 
of  the  configuration,  and  on  the  outer  boundary  of 
the  computational  domain,  before  constructing 
grids  in  the  interior  of  the  domain.  The  first 
step  is  to  convert  the  input  geometry  of  each 
component  of  the  configuration  to  a  continuous 
function  of  two  parametric  coordinates  using  the 
Coons  bl-cublc  patch  formulation.  This  trans¬ 
formation  from  a  discrete  to  a  continuous  descrip¬ 
tion  of  the  surface  allows  the  component  inter¬ 
sections  to  be  calculated  accurately.  The  surface 
grids,  with  a  topology  derived  from  the  topology 
defined  above,  are  then  generated  in  cerms  of  the 
parametric  coordinates  by  solving  a  3ec  of  non¬ 
linear  elllpcic  equations  as  proposed  by  Thompson, 
Thames  and  Mascln35,  using  over-relaxation  or 
approximate  factorisation  methods.  The  solution 
is  Chen  mapped  back  Co  physical  space  by  use  of 
the  bl-cublc  patches,  which  ensures  consistency 
with  che  original  georaecry,  and  provides  boundary 
conditions  for  the  generation  of  the  field  grid. 

Field  Grid  Generation 

The  approach  adopted  ac  ARA  Co  generate  the 
field  grid  has  been  to  solve  the  elliptic  Thompson 
equations  in  ch ree-d ime ns  Iona  1  form  on  a  block-by¬ 
block  basts  for  Che  entire  field.  The  physical 
locations  of  the  interior  block  boundaries  become 
fixed  in  che  course  of  an  iterative  solution.  The 
Jim  is  to  produce  a  grid  vlch  smooch  variation 


across  block  boundaries  as  well  as  within  a  block. 
The  presence  of  singular  polncs  and  lines  is 
accepced  (and  may  be  inevitable).  To  control  che 
distribution  of  grid  points  within  a  block  che 
control  functions  used  as  forcing  cerms  in  che 
Thompson  equations  are  expressed  in  terms  of  che 
point  distribution  on  a  given  boundary,  as 
suggested  by  Thomas  and  hiddlecof f 36 .  When  used 
with  mulciblock,  considerable  care  is  needed  with 
this  approach.  Grid  cross-over  can  occur  and 
point  distributions  frequently  fall  short  of  che 
quality  attainable  in  two-dimensional  grids.  Grid 
cross-over,  which  is  apparently  only  possible  with 
discretised  forms  of  the  Thompson  equations,  seems 
most  likely  to  occur  near  sharp  edges,  such  as 
wing  tips  and  cradling  edges. 

Automatic  Grid  Generation 

The  generation  of  a  satisfactory  grid  for  a 
practical  configuration  is  an  exceptionally 
demanding  Cask,  requiring  both  skill  and 
experience  and  absorbing  time  and  efforc.  It  is 
desirable,  therefore,  to  seek  means  of  automating 
the  task  as  far  as  possible,  while  maintaining  the 
quality  of  the  grid.  This  is  not  a  well-defined 
objective.  The  possible  level  of  automation,  what 
is  meant  by  a  satisfactory  quality  of  grid,  and 
the  nature  of  the  task  Itself  are  indefinite  and 
varying.  However,  some  guidelines  can  be 
identified. 

The  automation  of  topology  definition 
presencs  the  greatest  challenge.  This  is  the  area 
of  grid  generation  where  skill  and  experience  are 
of  tantamount  importance,  and  the  generation  of  a 
suitable  topology  is  a  necessary  condition  for  the 
generation  of  a  good  grid.  However,  topology 
definition  is  an  essentially  creative  task  -  for  a 
given  configuration  there  is  an  unlimited  number 
of  possible  topologies  -  and  so  automation  of  this 
process  is  particularly  difficult.  By  contrast, 
once  che  copology  has  been  defined,  the  generation 
of  che  surface  and  field  grids  is  relatively 
straightforward. 

Complete  automation  of  grid  generation  for 
arbitrary  configurations  does  not,  at  present, 
seem  feasible.  However,  for  particular  types  of 
configuration,  such  as  a  wing-body-canard,  the 
range  of  appropriate  topologies  can  be  limited  and 
it  becomes  feasible  to  develop  an  automatic 
routine  for  Copology  definition  within  this  subset 
of  topologies.  Nevertheless,  it  is  expected  that 
some  degree  of  user  interaction  will  still  be 
desirable. 

The  ARA  work28  on  grid  generation  for  wing- 
body-canard  configurations  illustrates  how  such  an 
approach  might  be  feasible.  Within  the  class  of 
wlng-body-canard  configurations,  there  Is  a  large 
range  of  'metrically'  different  configurations, 
associated  with  different  wing  sweeps,  different 
relative  location  of  canard  and  wing,  etc. 

Although  the  range  of  appropriate  topologies  is 
limited  there  are  nevertheless  many  possible  topo¬ 
logies  for  this  range  of  configurations-  In  order 
to  proceed,  ARA  choose  Co  have  C-grlds  around 
screamwlse  sections  of  lifting  surfaces,  an  O-grid 
around  the  body  section  and  an  H-grld  around  span- 
wise  sections  of  che  lifting  surfaces.  These 
grids  are  embedded  In  a  'global'  H-grld.  This 
provides  a  basis  for  automatic  copology  definition. 
The  user  must  supply  the  key  geometrical  features 
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chat  determine  the  topology,  such  as  whecher  the 
canard  is  to  lie  on  Che  same  coordinate  surface  as 
the  wing.  A  simple  format  is  provided  for  this  as 
a  'user  interface'.  In  addition,  the  user  speci¬ 
fies  the  number  of  blocks  in  each  of  the  coordi¬ 
nate  directions  of  the  global  H-grid.  The  number 
accually  used  may  thus  be  greacer  than  the  minimum 
required  for  a  given  topology,  for  example  to  give 
Increased  grid  density  in  appropriate  regions. 

The  surface  grids  are  generated  with  a  surface 
topology  derived  from  the  field  topology,  and 
require  two  inputs.  The  first  is  used  to  define 
the  grid  density  in  the  region  of  leading  and 
trailing  edges  of  the  lifting  surfaces,  while  the 
second,  from  the  topology  definition,  provides  for 
consistent  grid  densities  on  the  surface  and  on 
the  outer  boundary  of  the  computational  domain. 
Generation  of  Che  field  grid  is  Chen  a  straight¬ 
forward  process. 

Examples  of  Applications  of  CFD  Methods 
Analysis  of  Wlng-3ody-Canard  Design 

The  multiblock  grid  generation  scheme 
described  above  has  been  used,  in  conjunction  with 
the  BAe  method10  for  solving  the  Euler  equations, 
to  calculate  the  inviscid  compressible  flow  about 
a  wing-body-canard  configuration,  known  as  M165. 
This  is  a  research  model  that  has  been  the  subject 
of  a  number  of  wind-cunnel  tests. 

The  first  step  in  the  use  of  the  miltiblock 
scheme  is  the  definition  of  a  suitable  topology. 
Here  a  C-H-0  copology  is  adopted  with  an  O-grid 
around  the  body  section.  Many  other  reported  grid 
generation  techniques  for  wing-body  combinations 
use  H-grid  structures  or  C-grids  alone,  with  the 
result  chat  Che  grid  around  the  body  is  inadequate 
particularly  in  the  important  nose  region. 


leading  edge  are  poorly  represented,  indicating 
chat  the  grid  near  the  leading  edge  needs 
refining.  The  predicted  shocks  are  stronger  and 
further  aft  than  the  observed  ones.  However,  the 
incorporation  of  viscous  effects  would  produce  a 
weakening  of  che  predicted  shocks  and  forward 
movement  in  position  -  we  can  infer  therefore  chac 
the  presenc  inviscid  results  for  shock  strength 
and  position  are  consistent  with  che  measurements. 

The  second  comparison,  for  a  supersonic 
onset  flow  (Mach  number  1.2),  again  at  6  degrees 
incidence,  is  presented  in  Figure  5.  Here  che 
pressure  distribution  on  the  wing  shows  a  scrong 
trailing-edge  shock  chat  is  well-predicted  by  the 
method.  Again,  overall  agreement  becween  theory 
and  experiment  is  good. 

A  further  set  of  results  is  shown  in 
Figure  6.  Here,  the  method  has  been  used  to  pre¬ 
dict  che  effect  of  the  canard  on  che  wing  press¬ 
ures  at  M  «  0.9  and  6  degrees  incidence.  The 
expected  effect  of  the  canard  downwash  can  be 
seen  in  the  results,  namely  a  reduction  in  wing 
loading,  particularly  inboard  of  che  canard  tip 
position,  again  in  good  overall  agreement  with 
experimental  results. 

Design  of  an  Advanced  Civil  Fan  Rotor 

The  next  example  of  the  application  of  CFD 
methods  is  the  design  of  a  fan  rotor  for  a  civil 
turbofan  engine21.  The  method  used  comprises 
three  main  components: 

(1)  SI  -  an  inviscid  Euler  method  (based  on  the 
mechod  of  Denton) 

(li)  BYL2  -  an  integral  boundary  layer  method 


I 


Canard  configurations  can  present  appreci¬ 
able  problems  for  grid  generation  if  there  is 
substantial  spanwise  variation  in  the  distance 
between  che  canard  and  che  wing.  For  the  con¬ 
figuration  considered  here,  che  trailing  edge  of 
che  canard  and  che  wing  leading  edge  lie  in  the 
same  vertical  plane  ac  Che  body  side,  but  Che  tip 
of  the  canard  lies  well  upstream  of  the  wing  lead¬ 
ing  edge,  due  che  large  sweep  of  the  wing.  This 
spanwise  variation  of  relative  position  of  the 
canard  and  wing  Imposes  conflicting  requirements 
on  che  body-side  and  canard-tip  grids,  and  this 
dilemma  has  received  special  attention  In  the 
presenc  example  in  chac  two  disclncc  topologies 
for  che  region  of  che  canard  and  wing  have  been 
cried. 

Figure  3  shows  one  of  the  surface  grids 
used.  The  grid  has  49  points  around  Che  canard 
section,  with  8  spanwise  stacions,  while  the  wing 
was  represented  by  78  polncs  around  the  section 
and  15  spanwise  scaclons.  The  0-grld  structure 
around  the  body  Is  clearly  visible.  The  complete 
grid  comprised  some  88000  cells,  constituting  430 
b locks • 

Three  sets  of  results  are  now  presented.  The 
first  set  are  for  a  freescream  Mach  number  of  0.9 
and  Incidence  af  5  degrees.  Figure  4  shows  a 
comparison  of  predicted  and  measured  pressures  on 
the  wing  'the  tanard  was  not  pres  su  re-p  lo  t  ted  on 
the  model).  food  overall  agreement  has  been 
obtained,  although  che  sue t  Ion  peaks  near  che 


(Ill)  S2  -  a  streamline  curvature  throughflou 
analysis. 

These  three  elements  are  combined  in  an  Iterative 
manner  to  produce  a  quasl-three-dlmenslonal  tech¬ 
nique  for  analysing  fan-blade  performance. 

The  Euler  method  and  Integral  boundary  layer 
method  are  combined  (using  the  Inverse  coupling 
techniques  described  in  Reference  2)  to  give  a 
method  (referred  to  as  S1BYL2)  that  Is  capable  of 
■predicting  blade-to-blade  flows  Including  separ¬ 
ation  induced  by  shock  waves  or  marked  adverse 
pressure  gradients  (diffusion). 

The  blade-to-blade  calculation  takes  place 
on  a  specified  axisymmecric  stream  surface,  and 
incorporates  the  effects  of  rocatlon,  varying 
radius  and  variation  of  stream  tube  thickness  In 
the  axial  direction.  This  Information  Is  obtained 
from  the  throughflow  calculation  method  (S2), 
derived  assuming  axisymmetry.  A  number  of  stream¬ 
lines  along  a  radius  are  'tracked1  by  the  S2  cal¬ 
culation,  which  accounts  for  the  blockage,  momen¬ 
tum  loss  and  change  of  axial  momentum.  This  data 
is  obtained  from  the  blade-to-blade  calculation, 
and  so  a  complete  solution  Is  obtained  by  coupling 
che  S2  and  S18YL2  methods  in  an  iterative  manner. 

The  technique  has  been  used  as  a  'design  by 
analysis'  method  In  conjunction  with  a  family  of 
parametrically  defined  blade  profiles.  The 
suction  surface  of  the  blade  is  made  up  of  four 
circular  arcs,  whllsc  che  pressure  side  Is  con- 
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strucced  from  two  cubic  splines.  The  suction 
surface  shape  is  defined  from  aerodynamic  con¬ 
siderations,  and  the  lower  surface  shape  is  then 
chosen  primarily  to  give  a  structurally  feasible 
blade  section. 

The  transonic  fan  rotor  specified  for  the 
design  is  typical  of  a  modern,  high  bypass-ratio 
turbofan.  A  bypass  ratio  of  1.8  was  specified 
with  a  tip  speed  of  450  metres  per  second.  The 
inlet  Mach  number  varied  from  0.8  to  1.5  hub  to 
tip,  with  a  corresponding  variation  of  flow  angle 
of  35  Co  65  degrees. 

The  main  considerations  during  the  design 
process  were  to  avoid  choking  due  to  thickness  and 
close  pitch  ac  the  blade  rooc,  to  maximise  dif¬ 
fusion  without  flow  separation  at  the  mid-span, 
and  to  avoid  shock-induced  separation  at  the  tip. 
This  latter  requirement  meant  reducing  the  Mach 
number  ahead  of  the  shock  by  supersonic 
compression. 

Care  was  taken  to  ensure  that  the  three 
types  of  blade  section  chat  met  these  three  design 
considerations  could  be  'stacked'  to  produce  a 
realistic  blade  shape,  and  the  S2  calculation  was 
repeated  a  number  of  times  Co  ensure  the  sections 
were  matched.  The  S2  calculation  used  19  stream¬ 
lines,  while  8  radial  stations  were  chosen  for  the 
blade-co-blade  calculations- 

Typical  sections  from  the  resulting  design 
are  shown  in  Figure  7,  together  with  predicted 
flow  properties  at  the  design  condition.  Corre¬ 
sponding  results  for  an  earlier  conventional 
design  are  Included.  The  results  show  reduced 
pre-shock  Mach  numbers  for  the  new  design,  and 
boundary  layers  that  are  attached,  unlike  those 
for  Che  earlier  design.  This  leads  to  a  marked 
reduction  in  loss  of  total  pressure.  Overall,  the 
fan  was  predicted  as  being  some  27.  more  efficient 
than  previous  conventional  designs,  and  this  gain 


be  structured  and  thus  computational  efficient  but 
not  well  suited  for  some  configurations,  or 
unscruccured  and  chus  relatively  inefficient  but 
versatile,  or,  in  a  multiblock  framework,  partly 
structured  and  partly  unstructured.  For  strongly 
interacting  flows  the  use  of  vi scous-inviscid 
coupling  techniques  may  give  way  enclrely  to  solu¬ 
tion  of  the  Navler-Stokes  equations,  but  new 
couplings  may  be  deVised  that  successfully  exploit 
the  discincc  viscous-inviscid  scruccure  of  high 
Reynolds-number  flows  for  computational 
efficiency.  Finally,  large-scale  computing  may 
evolve  slowly,  with  reliance  on  the  present  super¬ 
computer  architecture  and  a  trend  towards  a  modest 
degree  of  multiprocessing,  or  be  transformed  by 
new  developments  in  the  simultaneous  use  of  large 
arrays  of  processors. 
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Fig  1  Mach  number  distributions  on  NACA  0012 
for  fitted  and  captured  shocks. 
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